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A  NON-LINEAR  SHOCK  WAVE  REFLECTION  THEORY 

ABSTRACT 

A  one-dimensional  theory  of  normal  reflection  of  blast  waves  from 
walls  is  given.  The  method  satisfies  the  initial  and  boundary  conditions 
of  the  problem.  It  is  shown  how  the  entire  reflected  wave  zone, 
including  the  reflected  shock  front  and  the  pressure  and  impulse  on  the 
wall,  can  be  calculated. 
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LIST  OF  SYMBOLS 


*  =  denotes  dimensional  quantities 

a  =  a  constraint  parameter  on  u 

b  =  a  constraint  parameter  on  u 

c*  -  sound  velocity  in  undisturbed  air 
oo 

c*  =  sound  velocity  in  region  traversed  by  incident  shock 
o 

e* 

o 

c  =  -5— 

O  C71 

OO 

c*  =  sound  velocity  in  region  traversed  by  reflected  shock 


C  =  ~— 

C* 

OO 

E*  =  specific  internal  energy  of  undisturbed  air 

E*  =  specific  internal  energy  in  region  traversed  by  incident  shock 


E 

o 


E* 


E* 

o 


E* 

00 


c*2 

00 


=  specific  internal 


energy  in  region  traversed  by  reflected  shock 


K*  _  2* 

00 


c* 

OO 


M*  =  mass  of  explosive 

p*^  =  total  pressure  of  undisturbed  air 

p*  =  total  pressure  in  region  traversed  by  incident  shock 


p*  =  total  pressure  in  region  traversed  by  reflected  shock 
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LIST  OF  SYMBOLS  (Contd) 


=  gas  constant 

=  specific  entropy  of  undisturbed  air 

=  specific  entropy  in  region  traversed  by  incident  shock 


S*  -  S* 
o  oo 


=  specific  entropy  in  region  traversed  by  reflected  shock 


S*  -  S* 


=  time 

c*  n* 

=  (  00*°°) 

=  particle  velocity  in  region  traversed  by  incident  shock 
=  particle  velocity  in  region  traversed  by  reflected  shock 


=  velocity  of  reflected  shock 


=  linear  distance 


P 

_  /  o° 

V  2  • 

c*  M* 
oo 


=  specific  heat  ratio  of  air 


,Y  - 

u  = 
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LIST  OF  SYMBOLS  (Conta) 


P*  =  density  of  undisturbed  air 
oo 

p*  =  density  in  region  traversed  by  incident  shock 


p*  =  density  in  region  traversed  by  reflected  shock 


YO  * 
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INTRODUCTION 


The  study  of  damage  to  structures  by  blast  waves  requires  analysis 
of  the  normally  reflected  pressure  on  the  wall.  To  solve  this  problem 
the  hydrodynamical  equations  of  flow  must  be  solved  in  the  region  between 
the  wall  and  the  reflected  shock  front.  Since  the  system  of  non-linear 
partial  differential  equations  of  flow  are  presently  solvable  only  by 
numerical  methods  that,  when  reliable,  are  somewhat  cumbersome,  and  the 
possibility  of  exact  analytical  solution  is  remote,  a  method  that  simpli¬ 
fies  the  mathematics  somewhat  and  is  more  amenable  to  analytic  solution 
is  desirable.  Here,  we  examine  a  direction  of'  simplification  that  reduces 
the  system  of  equations  for  one -dimensional  flow  to  a  system  of  ordinary 

differential  equations,  by  introduction  of  a  constraint  that  satisfies 

1* 

the  initial  and  wall  conditions,  as  discussed  oy  Makino. 

2 

Makino  and  Shear  have  obtained  a  theory  of  reflected  impulse  at  the 

wall  by  regarding  each  element  of  the  incident  wave  to  be  individually 

reflected  like  the  shock  front,  but  this  is  a  zero'th  order  approximation, 

3 

not  satisfying  derivative  conditions  at  the  wall.  Chang  and  Laporte 
have  obtained  two  theories  of  shock  reflection.  One  is  series  expansion 
about  the  point  of  reflection,  which,  if  truncated  for  practical  pur¬ 
poses,  may  not  fully  satisfy  the  wall  conditions.  The  other  theory,  which 
assumes  the  particle  velocity  to  be  zero  all  along  the  reflected  shock 
line,  may  also  not  fully  satisfy  the  wall  conditions.  Also,  both  theories 
are  for  the  calculation  of  the  reflected  shock  line  only.  In  the  theory 
we  present  here,  we  consider  calculation  of  the  entire  reflected  wave  zone, 
from  the  shock  front  through  interior  points  to  points  on  the  wall,  such 
that  wall  conditions  are  satisfied  through  certain  derivatives.  However, 
it  is  probably  most  useful  only  in  that  phase  of  the  wave  exerting  the 
greatest  stress  on  the  wall,  which  is  the  part  of  greatest  interest  for 
damage  studies . 


Superscript  numbers  denote  references  which  may  be  found  on  page  26. 
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Ryzhov  and  Khri stianovich  have  developed  a  theory  on  the  problem 

of  two  dimensional  regular  reflection,  but  the  theory  assumes  isenlropic 

flow  and  is  therefore  applicable  to  weak  shocks  only,  and  further,  the 

.  .  ,  . 5  , 

boundary  conditions  are  approximated,  smnctiapin  nas  improved  tne 
theory  with  respect  to  boundary  conditions,  but  has  not  extended  the 
theory  to  non-isentropic  flow. 

C  r7 
\J  j 

While  the  so-called  self-similar  type  solutions  can  be  extended 
to  the  reflection  problem  as  an  approximation,  the  choice  of  the 
similarity  form  to  be  assumed  is  made  difficult  by  the  strong  influence 
of  the  wall  and  by  the  non-constancy  of  the  quantities  in  front  of  the 
reflected  shock. 

The  example  considered  here  is  for  plane  flow,  or  spherical  flow 
at  distances  sufficiently  far  from  the  center  of  energy  release  that 
planar  approximation  suffices.  It  is  shown  how  the  flow  parameters 
behind  the  reflected  wave,  in  particular  the  pressure  on  the  wall  as 
a  function  of  time  and  the  impulse,  can  be  obtained.  For  cylindrical 
and  spherical  flows,  the  same  method  is  applicable  if  the  reflecting 
surfaces  have  the  corresponding  symmetries. 


FLOW  EQUATIONS 


The  non-dimensionalized  equations  of 

rS  t  mo  n  c  l  on  a  1  mo  +  -i  on  of1  a  1  r»  a  no  *1  n  "F/n  “I  o  n'T  an 

VA.  UiV  *10  -1-  VUU-L  U  x  V*i  i-  IU.-4.  A  j  -kli  iJUXWl  J.  WA,  i.  J. 

of  mass 


flow  describing  the 
..  .  8 

ooo  n  H  i  n  T  o  c  oonco 

v-/  W  i.  UXiiUV1-  O  3  ^  ij  W 


one- 


m  r  -i  on 

A  V  Wii 


( la) 


n.  +un  +ocu  =0. 

‘t  x  X  X 


conservation  of  momentum 


(lb) 

adiabaticity 

(lc) 


7  p^  +  u+  +  uu  =  0  , 
p  X  o  X 


S  +  uS  =  0 , 
t  x 
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where  t  is  the  time  with  the  non-dimensionalizing  scaling  factor  (ambient 
sound  speed  times  ambient  pressure /mass  of  explosive)  ,  x  is  the  distance 
with  the  non-dimensionalizing  scaling  factor  (ambient  pressure/( ambient 
sound  speed)  /mass  of  explosive)  ,  p  is  the  pressure  in  units  of  the 
ambient  pressure,  u  is  the  particle  velocity  in  units  of  the  ambient 
sound  speed,  c  is  the  sound  velocity  in  units  of  the  ambient  sound  speed, 
p  is  the  specific  heat  ratio  Y  (assumed  constant)  times  the  density  in 
units  of  the  ambient  density,  and  S  is  the  excess  entropy  over  the  ambient 
in  units  of  the  gas  constant  R. 

This  system  of  equations  is  supplemented  by  the  equation  of  state, 
which,  for  illustrative  purpose,  we  assume  to  be  ideal: 


(2a) 

1/Y  t  Y  -  1  o \ 

P  -  YP  exp  (-  —  S)  , 

(2b) 

E  1  P .  1  .. 

Y  -  1  P  Y(Y  -  1)  ' 

where  E  is  the  non 

-dimens i on ali zed  energy. 

The  Eulerian  coordinates  x,  t  in  the  equations  above  are  replaced 

by  Lagrange  coordinates.  We  define  m  to  be  the  mass  integral 

(3) 

m  =  j"  p(x,t)  dx  , 

where  the  integration  is  performed  on  a  constant  t  line  starting  from  the 
wall.  From  this  definition  and  from  the  continuity  Equation  (la)  in  the 
form 


(M 


°t  '  -(pu)x  • 


we  obtain 

( 5a) 

(5b) 


mx  =  P  ’ 
mt  =  -pu  , 
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and  also 


(6)  x+(m,t)  =  u. 

G 

Using  Equations  (5)  and  (2),  we  put  Equation  (l)  in  the  form 


(7) 


1  0  0  0  0  0 

0  10  ypp  0  0 

0  0  1  0  10 


V  = 


where  V  s  col  j|  St ,  pt ,  ut ,  V  pm, 
differentiation  of  Equation  (7)  with  respect  to  rn  and  t  gives 
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The  Rankine-Hugoniot  conditions  that  must  be  satisfied  across  the 
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shock  front  are,  in  dimensionless  form, 
conservation  of  mass, 


(10a) 


p(U  -  u)  =  p  (U  -  u  )  , 
o  o 


Ik 


conservation  of  momentum. 


(,10b; 


_  /  T  Y 

p^u 


-  u)‘ 


=  P  (u  - 
o 


n2 

u  ) 
o 


conservation  of  energy, 


(10c) 


hu  -  u)2  +  E  +  £  =  i(U  -  u  )L  +  E  +  — 

d  p  d  O  O  p 


d.x 

where  U  is  the  shock  velocity  —7  . 

These  conditions  simplified  by  Equation  (2)  give 


(Ha) 


u  =  u(p;uo,po>So)  =  u_ 


+  (P  -  P0} 


exp 


(^s  ) 

Y  O 


,(Y  +  1}  pi/Y(P  +  u2  P  i 


1/2 


(lib) 


S  =  S(p;u  p  S  )  =  S 

OOO  o 


V  .  ,  P,i/T  /  U"P  +  Po  \ 
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(  1  1 
\  -L-l.'-  / 


u  =  U(p;uo,po,So)  -  uo 


Y  -  1 


(p  +  V  P  )  exp  ( 

- 2 _ I 

/,  2,  1/y 

(1  -  p  )  p  1 


S_) 


1/2 
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3  3 

Let  D  =  —  +  n  l  II  -  u)  —  be  the  differential  ODerator  in  the  direct  ion 

at  " ‘  am  1 

of  the  shock  line.  Implicit  differentiation  of  u  and  S  along  the 
reflected  shock  line  gives 


C  i  P  a  ) 


3U 

3p 


_  /  T  I  ,  ,  \  , 

p  v  u  —  u/u 


m 


_  /  ;  i  ■■ )  r 

P^U  "  u  dp"  ^m 


dU 


DP, 


dU 


Du  + 


dU 


DS 


(12b) 


n  —  n  ( F  [  —  1 1  i  n  +  n  ( U  —  u)S  =F 

3p  *”t  “/  3p  ■  PV“  1  m  2  ’ 


f-Dp  +  |LDu  +||_DS 

3p  o  3u  o  3S  o 


where  the  shock  path  is  given  by 


(13) 


Dm  =  p ( U  -  u)  . 


The  quantities  ,  identified  by  subscript  o,  that  result  from  the  incident 
wave  we  assume  to  be  known  as  functions  of  m  and  t.  From  Equation  ( 11 ) , 
the  first  partial  derivatives  with  respect  to  the  arguments  p,  uq ,  p^ , 

S  ,  become 


fii,h 

v  ita  / 


3u 
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/ 
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- 
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PC 
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(14b) 
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=  i  , 


(14C) 
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1 6 
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CONSTRAINTS 

For  simplicity,  it  is  assumed  that  the  reflecting  wall  is  at  m  =  x  =  0 
(see  Figure).  At  the  wall  the  particle  velocity  must  satisfy  the 
conditions 


(15) 


/  anu(x,t) 

l 


at 


) 


/  x  =  0 


From  this  result  and  the  ideal  gas  law  (2a),  we  can  show  by  taking  higher 
derivatives  of  (la)  and  (lb)  that  u  must  satisfy  also  the  condition 


(1 6) 


( a2u( x,t ) ^ 

V  ax2  /  x  = 


n 


V 


0  s 


IT 


We  now  impose  some  constraints  on  the  flow.  There  is  an  infinity 
of  ways  of  doing  this.-  The  choice  will  depend  on  the  nature  of  the 
ambient  conditions  and  the  incident  wave.  For  blast  waves  in  still  air, 
in  the  neighborhood  of  the  wall  at  all  times  and  also  in  the  neighborhood 
of  the  asymptotic  shock,  we  expect  the  particle  velocity  to  vary  slowly 
with  x,  and  so  we  choose 

(17)  u  =  a(t)x  +  b(t)x  , 


where  a(t)  and  b(t)  are  functions  of  t  only.  This  expression  satisfies 
both  boundary  conditions  (15)  and  ( l6 )  ,  and  also  permits  the  initial  and 
shock  conditions  to  be  satisfied. 


For  subsequent  purpose,  we  differentiate  Equation  (17)  twice  with 

“h  /~\  m  « 

V—  VJ  u  v  Ail  • 


(18a) 


u  = 


u  + 


(18b) 
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u  =  ^  -  —  [p  -  (y  —  1)  PS  ] 
mm  2  yp  m  m 
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DETERMINATION  OF  REFLECTION  DOMAIN 

Equations  (l8a),  (12),  and  (7)  are  six  equations  involving  the  nix 
components  of  V.  Solving  for  V,  we  obtain 


(19) 


where 


V  = 


-  H{P ( U  -  u)t~  +  1} 
3p 


-G  -H 


—  (p(U  -  u)f^  +  1} 

YPP  3p 


r. i  -l. 

'  Ur. 
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P(U  -  u) 


[f2  *  §  "*<“  - 


U)  - 


H} 
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\  v  —  U J  n.  ’  -L  J 

3p 


H  = 


Yp(u  +  2bx~J 

x{p  (,  U  -  uj—  +  1} 
3n 


By  differentiating  Equation  (12)  along  the  reflected  shock  line,  we  have 


(20a) 
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C  20  b ) 
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Equations  (20),  (l8b),  and  (8)  are  nine  equations  involving  the 
nine  components  of  W.  Solving  for  W  gives 


(22) 


w  =  b_1  y  , 
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Solving  Equation  ( 18a  )  for  b(t),  there  is  obtained 


(2h) 


b(t)  = 


px  u  -u 
m 
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Differentiating  this  expression  with  respect  to  t  along  the  reflected 


(25) 
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Equations  (2a),  (ll),  (13),  (25),  together  with  the  implicit  derivative 

(26) 

give  the  system  of  ordinary  differential  equations 
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Using  (19) 

and 

(22)  , 

we 

can 

solve 

this 

system. 

The  solutions  of  this  system  give  b(t)  as  a  function  of  time;  a(t) 
is  then  determined  by 


I  \ 
;  cuj 


-  _  u  -  b(t )x“ 

av  uj  -  - 
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The  flow  behind  the  reflected  shock  can  now  be  calculated  as 
solutions  of  ordinary  differential  equations,  with  each  point  on  the 
calculated  shock  line  serving  as  an  initial  point,  x  is  given  by 


(29) 


x  =  u  =  a(t)x  +  b(t)x~ 


which  integrates  along  the  m  =  constant  line  into 


(30) 


x 


t(m) 


a(t  ^  at  \ 


_  C  ^  2  x2  ( ir. )  b  ( t ) 

-  t ( m) 


exp 


2a(i)dTldt  +  1 


-1/2 


where  x(m) ,  t(m)  is  a  point  on  the  shock  line.  This  value  of  x 
substituted  into  Equation  (29)  gives  u  as  a  function  of  m,  t.  The 
pressure  is  then  obtained  from  Equation  (7): 


(31) 


Pt  =  -  YPPum  = 


=  -  YP [a(t )  +  3b (t )xc ] 


which  integrates  into 


(32) 


p  =  p(m)  exp  Ut  Y  [a(t )  +  3b(t)x2]dti  . 


t(m) 


On  this  m  =  constant  line,  the  entropy  S  =  S(m)  is  constant  and  is  given 
on  the  shock  line.  The  density  then  follows  from  Equation  (2). 

On  the  wall  the  pressure  simplifies  to 


(33) 


p(o,t) 


p(0) 


l  /  ri  \ 

rvu; 


2k 


The  positive  impulse  imparted  to  the  wail  we  take  to  be 


/  oi.  \ 
v  _>*+  / 


Pos .  impulse 


t{0) 


p(0)  exp  /-  f  ^  a(t)dt] 


1  J 


t  ( 0 ) 


where  t  is  the  time  at  which  p(0,t)  =  1. 


RALPH  E.  SHEAR 


RAY  C. 


-  1 1  dt 
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